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A coherent description of the anomalous Hall effect (AHE) is presented that is applicable to pure 
as well as disordered alloy systems by treating all sources of the AHE on equal footing. This is 
achieved by an implementation of the Kubo-Stfeda equation using the fully relativistic Korringa- 
Kohn-Rostoker (KKR) Green's function method in combination with the Coherent Potential Ap- 
proximation (CPA) alloy theory. Applications to the pure elemental ferromagnets bcc-Fe and fcc-Ni 
led to results in full accordance with previous work. For the alloy systems fcc-Fe^Pdi-a; and fcc- 
Ni^Pdi-a; very satisfying agreement with experiment could be achieved for the anomalous Hall 
conductivity (AHC) over the whole range of concentration. To interpret these results an extension 
of the definition for the intrinsic AHC is suggested. Plotting the corresponding extrinsic AHC versus 
the longitudinal conductivity a linear relation is found in the dilute regimes, that allows a detailed 
discussion of the role of the skew and side-jump scattering processes. 
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During the last years the anomalous Hall effect (AHE) 
has received great interest. This is partly caused by 
its close connection to the spin Hall effect (SHE), that 
possesses a large potential for application in the rapidly 
growing field of spintronics [l|. On the other hand, many 
theoretical investigations are devoted to the development 
of a coherent description of these quite complex phenom- 
ena 0. 

As was already pointed out by Karplus and Luttinger 
3 the ultimate origin for the AHE in ferromagnets is the 
spin-orbit coupling (SOC) that - together with the spon- 
taneous magnetization - leads to a symmetry breaking. 
As was demonstrated by experiment 0, Q and is obvi- 
ous from the work of Karplus and Luttinger the AHE 
is present even in pure systems. This so-called intrinsic 
AHE could later be connected to the Berry-phase Q and 
corresponding ab initio results could be obtained during 
the last years using an expression for the anomalous Hall 
conductivity (AHC) cr xy in terms of the Berry curvature 
For diluted and concentrated alloys, on the other 
hand, the occurrence of the AHE was primarily ascribed 
to the spin-dependent skew or Mott [3, [l(| and the so- 
called side-jump [ll| scattering mechanisms. The latter 
one is caused by the anomalous velocity, a first-order rela- 
tivistic correction to the non-relativistic velocity operator 
connected to SOC. Interestingly, scaling laws connecting 
the AHC <7 xy and the longitudinal conductivity er xx (see 
below) could be derived for these two extrinsic mecha- 
nisms 0] ■ Their treatment in connection with a descrip- 
tion of electronic transport in terms of wave packet dy- 
namics was discussed in detail recently by Sinitsyn [12j . 
When dealing with the extrinsic AHE in disordered sys- 
tems, however, disorder was treated so far o nly by model 
potentials [13[ or by a damping parameter [l5[ . 
Crepieux and Bruno [l6j performed qualitative investi- 



gations on the AHE on the basis of the Kubo-Stfeda 
equation. This equation is derived from Kubo's linear 
response formalism supplying a suitable basis for investi- 
gations based on a realistic description of the underlying 
electronic structure (see Ref. [l7l and below). An alterna- 
tive description of the AHE with a wider regime of appli- 
cability is achieved by using the non-equilibrium Green's 
function formalism. Using a suitable, but still tractable, 
model description for the electronic structure Onoda et 
ai. plEl could divide the range of c xx covered typically 
by real materials into three regimes with different scaling 
laws connecting cr xy and <r xx . 

In this communication results for the AHC obtained 
using the Kubo-Stfeda equation are presented. Using a 
fully relativistic Green's function formulation in combi- 
nation with a reliable alloy theory a coherent description 
for pure as well as diluted and concentrated alloys could 
be achieved that treats intrinsic and extrinsic sources of 
the AHE on equal footing. 

The Kubo linear response formalism supplies an ap- 
propriate basis to deal with electronic transport in mag- 
netic metallic systems. Making use of a single-particle 
description of the electronic structure and restricting to 
the case T = K one is led to the Kubo-Stfeda equa- 
tion for the electrical conductivity tensor er [l8| . For cu- 
bic systems with the magnetization along the z-direction, 
the AHE is described by the corresponding off-diagonal 
tensor element or anomalous Hall conductivity er xy given 
by: pi [H] 



Trace (j x (G+ - G" %G~ 
-] x G+j y (G+ -G~)) c 
Trace ((G+ - G")(r x j y - f y j x )} c (1) 
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Here Q is the volume of the unit cell, N is the num- 
ber of sites, while f and j are the position and current 
density operators, respectively. For the cubic systems 
considered here the last term is site-diagonal for symme- 
try reasons. As furthermore all system considered here 
are metallic it has been omitted [I?} • The electronic 
structure of the system is represented in terms of the 
single-particle retarded (G + ) and advanced (G~ ) Green's 
functions at the Fermi energy Ep. Within the present 
work these functions have been evaluated by means of 
the multiple scattering Korringa-Kohn-Rostoker (KKR) 
formalism [2(3] • The chemical disorder in the investigated 
random substitutional alloys has been accounted fo r by 
using the Coherent Potential Approximation (CPA) |21| . 
This alloy theory supplies a reliable framework to per- 
form the configurational average indicated by the brack- 
ets (...} c in Eq. (JT|). It includes, in particular, a clear 
definition for differences of configurational averages like 
(jxG + jyG~) c — (jxG + ) c (jyG~)c- These so-called vertex- 
corrections correspond to the scattering-in terms within 
semi-classical Boltzmann transport theory [22| . 

Dealing with the AHE requires to account for the 
influence of spin-orbit coupling in an appropriate way. 
This is achieved by using the four-component Dirac 
formalism. 23] In combination with spin-density func- 
tional theory in its local approximation (LSDA) the cor- 
responding Dirac Hamiltonian is given by: |24j 

H D = col ■ p + fimc 2 + V + (3Z Z B . (2) 

Here p = — i?iV is the canonical momentum operator, a 
and j3 are the standard Dirac matrices, [23[ while V and 
B represent the spin-independent and spin-dependent, 
respectively, effective LSDA potentials for the magne- 
tization along z. Within the fully relativistic frame- 
work adopted here the current density operator j is given 
by:[H 

j = -c|e|a. (3) 

To allow for a more detailed discussion on the origin of 
the AHE it is useful to introduce the alternative current 
density operator: [25j 

Jp = rF7r(p+ — a + -PZ z (a x ,a v ,0) T J ,(4) 

TO + E/C 1 \ c c ) 

that is equivalent to j given by Eq. @ and that can 
be derived from the anti-commutator of j and the Dirac 
Hamiltonian T-Lu given by Eq. ([2]). 

Recently, the intrinsic AHE of the pure ferromagnets 
F^Co and Ni @,H,[lil as well as ordered FePt and FePd 
[271 ] has been investigated theoretically on an ab initio 
level using the formulation for t7 xy in terms of the Berry 
curvature. Alternatively, the tensor element <r xy can be 
obtained directly from the expression given in Eq. ([1]) 
that is evaluated by Fourier transformation leading to a 
corresponding Brillouin zone integration 22] . As the in- 
tegrand shows a <5-function like behavior for pure systems 



TABLE I: The intrinsic AHC of bcc-Fe and fcc-Ni from ab 



initio theoretical as well as experimental (Exp.) investiga- 
tions. 



<r xy (mil cm) 1 


bcc Fe 


fee Ni 


present work 


0.638 


-1.635 


Yao et al. [7] 


0.753 




Wang et al. [8] 


0.751 


-2.203 


Yao [28] 




-2.073 


Exp. [4, 5] 


1.032 


-0.646 



a small imaginary part e has to be added to the Fermi 
energy Ep and an extrapolation to zero has to be made 
for e. For the calculations of <7 xy performed for bcc-Fe 
and fcc-Ni e has been varied between 10 -3 and 10~ 6 Ry. 
To ensure convergence of the Brillouin zone integration 
about 10 9 k-points have been used. The resulting AHC 
of bcc-Fe and fcc-Ni is given in Table I together with 
experimental data as well as results of previous ab initio 
work. fA \§L I2H Taking into account that SOC was treated 
by approximate schemes within the corresponding calcu- 
lations, the agreement is quite satisfying. 

The expression for er xy in terms of the Berry curvature 
used within previous work is completely equivalent to the 
Kubo-Stfeda equation used here, as both approaches are 
based on Kubo's linear response formalism and adopt a 
single-particle description for the electronic structure for 
T = K. 17] However, it should be stressed that the 
Berry curvature is usually formulated in terms of Bloch 
states implying translational symmetry this way. Calcu- 
lations of er xy for disordered alloys with broken transla- 
tional symmetry are therefore not possible on this ba- 
sis while the Kubo-Stfeda equation supplies an adequate 
framework for such investigations. As the k-dependent 
integrand connected with Eq. ([T|) gets smeared out in this 
case due to the chemical disorder it looses its ^-function 
like behavior. For that reason, broadening via a complex 
Fermi energy with e > is not necessary for calculations 
on alloys. 

Within the present work, corresponding calculations 
have been done for the alloy systems fcc-Fe x Pdi_a; and 
fcc-NLrPdi^a;. Both systems are formed by the nearly 
ferromagnetic transition metal Pd with an elemental 3d- 
ferromagnet leading to a very low critical concentration 
x cr it for the on-set of spontaneous ferromagnetic order 
(x^ [t ra for Fe.Pdi^ and x™ t « 0.02 for Ni s Pdi_ x ). 
As can be seen in Fig. [1] calculation of the AHC via 
Eq. ((T|) leads to a very satisfying agreement with the 
experimental data [29j . in particular, for Fe^Pdi-j;, that 
are available over a wide range of composition. In partic- 
ular the change in sign of cr xy with composition observed 
for both alloy systems is well reproduced by the calcula- 
tions. As one might speculate from the data for elemental 
bcc-Fe and fcc-Ni in Table I the sign of c xy on the Pd- 
poor side of fcc-Fe^Pdi-a is indeed positive while it is 
negative for fcc-Ni a .Pd 1 _ a; . For the Pd-rich the situation 
is reversed, clearly showing that the elemental ferromag- 
net is the primary source for the AHE in these two alloy 
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FIG. 1: (Color online) The AHC of fcc-Fe^Pdi-* and fcc- 
Ni^Pdi^a;. The total AHC o" xy (full squares) has been cal- 
culated including the vertex corrections (a x 
the intrinsic AHC a 
omitting them (a\ 



ntr 
xy 

ntr _ 

xy ~~ 



cr^p), while 



(open squares) has been obtained by 
o"xy VC ). In addition, experimental 
data [29j ] for er xy (full circles) determined at T = 4.2 K are 
shown. Further, an estimation for the side-jump contribution 

S J tn the PYfrinsir AflT nf frv-Fp _PH, _ and fcC-Ni^Pdl-^ 

is shown (open 

diamonds; see text). 



erjy to the extrinsic AHC of fcc-Fe^Pdi 
calculated as the difference cr X y tr — <r X y tr [V 



xy 



systems (see below). 

To get a more detailed insight into the mechanism re- 
sponsible for the AHE in the investigated alloys, a decom- 
position of the AHC has been performed. A formal basis 
for this is provided by the representation of the Kubo- 
Stfeda equation in terms of Feynman diagrams. jl6| From 
this it can be seen that the skew and side-jump mecha- 
nisms are exclusively connected to diagrams involving the 
vertex corrections. The remaining diagrams are standing 
for products of the type (jxG + ) c (j y G~) c that correspond 
to the intrinsic AHE and correction terms due to chem- 
ical disorder. It seems therefore sensible to extend the 
definition of the intrinsic AHC <7™ tr to the case of di- 
luted and concentrated alloys by combining all contri- 
butions not connected to the vertex corrections. This 
obviously allows to calculate the total AHC (cr xy ) and 
the intrinsic one (cr X y tr ) by evaluating the Kubo-Stfeda 



equation (Eq. (JTJ) with and without, resp., including 
the vertex corrections (VC); i.e. identifying cr xy 



and cr X y tr = o- x ° vc , respectively. As seen in Fig. [TJ cx™ tr 
gives a major contribution to the total AHC er X y of fcc- 
Fe x Pdi_ x and shows in particular also a change in sign 
with varying concentration. For fcc-NLcPdi-^, on the 
other hand, CT xy tr varies weakly with composition and ex- 
trapolates rather well to the intrinsic AHC of pure Ni 
(see Table I). For both alloy systems cr x " tr ss 1 (mficm) -1 
when xpd approaches 1 indicating that the intrinsic AHC 
is primarily determined by the properties of the Pd-host 
in the dilute regime. These findings obviously justify the 
extension of the definition for cr xy tr to represent all con- 
tributions not connected to the vertex corrections. 

The longitudinal conductivity <7 XX of fcc-Fe^Pdi-j; 
and fcc-Ni^Pdi-z lies nearly exclusively in the so-called 
super-clean regime with cr xx > (/ificm) -1 .[l4],[l5| For this 
regime the skew scattering mechanism should dominate 
<7 xy obeying the relation cr xy = S er xx , with S being the so- 
called skewness factor [3, [l5[ . Accounting for all three 
mechanisms one is therefore led to the decomposition: Q 



'xy 



-xy 



S a x 



S3 = ^mtr 

^ xy xy 



xy 



(•5) 



that may be seen as a definition for the side-jump con- 
tribution cr xy .0] In fact, a plot of cr xy versus er xx with 
the concentration as an implicit parameter was used in 
the past to decompose the experimental AHC of alloy 
systems accordingly. [3GH32] ] 

In Fig. [5] the extrinsic AHC of fcc-Fe^Pdi-a; and fcc- 
Ni^Pdi-j; defined as cr X y tr = er xy — cr xy tr is plotted versus 
the longitudinal conductivity <t xx . Obviously, the rela- 
tion suggested by Eq. (J5]) is well fulfilled on the Pd-rich 
side of both systems as well as on the Ni-rich side of 
Ni^Pdi-a;. Extrapolating for these regimes to er xx = 
allows to deduce the corresponding skewness parame- 
ters and side-jump term cr x J y (Fe^Pdi-^ for xpd > 0.9: 

w 0.1 (mOcm)" 1 ; Ni^Pdi-z for 
3 and fJ x J y ~ 0.4 (mficm) -1 and 
" 3 and oil ~ A Q '™rw,,-:n 



and cr x J y 
2.0 • 10- 



5 = -2.7- 10 
x Pd > 0.9: S 

for x m > 0.9: S = -6.6 • lO" 3 and a^ y « 4.8 (mSlcm)- 1 ) 
These results show clearly that the skew scattering mech- 
anisms by far dominates er xy tr in the dilute regimes. For 
the two alloy systems the corresponding skewness fac- 
tor S is found comparable in magnitude but different in 
sign on the Pd-rich side (see above). This once more 
demonstrates that the skew scattering mechanism has to 
be associated primarily with the solute component Fe or 
Ni, respectively. 

As emphasized above, Eq. ([5]) can be seen as a def- 
inition for various extrinsic contributions to a. 



xy 



ac- 



cording to their scaling behavior. An alternative way 
to define the side-jump term er x J y is to make use of its 
connection with the anomalous velocity, that is a cor- 
rection to the non-relativistic current density operator 
j nr = — T^ffV. Within the relativistic approach used 
here, an estimate for cr x J y can be made using the alterna- 
tive current density operator j p with the potential terms 
V and B suppressed (see Eq. The corresponding ex- 
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FIG. 2: (Color online) The extrinsic AHC <r X y tr versus a xx for 
fcc-Fe^Pdi-a; and fcc-NLjPdi-x. The straight lines represent 
extrapolations of the data for xp^ > 0.9 (Pd-rich) and x-si > 
0.9 (Ni-rich) to cr xx = 0. 



trinsic AHC <7^ tr [V] = cr xy [V] - 4 n y tr [V] allows to write 
a% « a%[V] = - <f'[V]. The results for a% ob- 



tained this way for fcc-Fe^Pdi-z and fcc-Ni^Pdi-^ are 
also shown in Fig. [TJ As one notes there is obviously a 
non-negligible concentration dependency for both alloy 
systems in particular on the Pd-rich side. In both cases, 
however, the numerical results are much smaller than for 
the intrinsic as well as the skew-scattering contributions. 
While this once more supports the conclusion that the 
AHE of the investigated alloy systems is dominated by 
the latter mechanisms, it also shows that the quantitative 
results for the side-jump term <7 x J y may depend strongly 
on the definition used. 

In summary, a coherent description of the AHE for 
pure metals and diluted as well as concentrated alloys 
on an ab initio level was presented based on a fully rela- 
tivistic implementation of the Kubo-Stfeda equation us- 
ing the multiple-scattering or KKR formalism in combi- 
nation with the CPA alloy theory. The intrinsic AHC 
obtained this way for bcc-Fe and fcc-Ni was found in sat- 
isfying agreement with previous ab initio work using an 
equivalent expression for cr xy in terms of the Berry cur- 
vature. Corresponding calculations for the alloy systems 
fcc-FezPdi-a; and fcc-Ni^Pdi-a, reproduced the available 
experimental data very well. Identifying the contribu- 
tions to er xy that are not connected to the vertex cor- 
rections with the intrinsic AHE of an alloy allowed to 
decompose the remaining extrinsic AHE. Plotting cr xy 
versus cr xx it was found that the skew scattering term 
by far dominates the side-jump contribution in the di- 
lute alloy regime. This conclusion could be supported 
by model calculations that supplied an estimate for the 
contribution to er xy due to the anomalous velocity. 
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